Charmonium properties in deconfinement phase in anisotropic lattice QCD 
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J/ty and n c above the QCD critical temperature T c axe studied in anisotropic quenched lattice 
QCD, considering whether the cc systems above T c are spatially compact (quasi-)bound states or 
scattering states. We adopt the standard Wilson gauge action and 0(a)-improved Wilson quark 
action with renormalized anisotropy a s /ot = 4.0 at j3 = 6.10 on 16 3 x (14 — 26) lattices, which corre- 
spond to the spatial lattice volume V = L 3 ~ (1.55fm) 3 and temperatures T ~ (1.11 — 2.07)T C . We 
investigate the cc system above T c from the temporal correlators with spatially-extended operators, 
where the overlap with the ground state is enhanced. To clarify whether compact charmonia survive 
in the deconfinement phase, we investigate spatial boundary-condition dependence of the energy of cc 
systems above T c . In fact, for low-lying S-wave cc scattering states, it is expected that there appears 
a significant energy difference AE = £(APBC) - £(PBC) ~ 2^/m 2 + 3tv 2 /L 2 - 2m c (rn c : charm 
quark mass) between periodic and anti-periodic boundary conditions on the finite-volume lattice. 
In contrast, for compact charmonia, there is no significant energy difference between periodic and 
anti-periodic boundary conditions. As a lattice QCD result, almost no spatial boundary-condition 
dependence is observed for the energy of the cc system in J/'fy and n c channels for T ~ (1.11 — 2.07)T C . 
This fact indicates that J/\t and r\ c would survive as spatially compact cc (quasi-)bound states be- 
low 2T C . We also investigate a P-wave channel at high temperature with maximally entropy method 
(MEM) and find no low-lying peak structure corresponding to Xcl at 1.62T C . 

PACS numbers: 12.38.Gc, 12.38.Mh, 14.40.Gx, 25.75.Nq 



I. INTRODUCTION 

To complete the phase diagram of quantum chromo- 
dynamics (QCD) is one of the most challenging attempts 
in particle physics. The difficulty of QCD originates 
from the nonperturbative nature in the low-energy re- 
gion, where the running coupling constant becomes large. 
As a consequence, color confinement and chiral symmetry 
breaking occur as nonperturbative phenomena, and the 
vacuum becomes the hadronic phase. On the other hand, 
at high temperature or high density region, color decon- 
finement and chiral symmetry restoration are expected to 
be realized. This phase is called the quark-gluon-plasma 
(QGP) phase. Actually, lattice QCD simulations show 
color deconfinement 1] and chiral symmetry restoration 
Q above the QCD critical temperature T c . The QGP 
phase transition is also investigated in various effective 
models 0, S IE IE S III 

As an important signal of QGP creation, J/ 1 ^ suppres- 
sion 0, [Hi LLjJ |l4jl3| was theoretically proposed in the 
middle of 80 's (HQ- The basic assumption of J/^> sup- 
pression is that J /\& disappears above T c due to vanishing 
of the confinement potential and appearance of the De- 
bye screening effect, which are actually shown in lattice 
QCD simulations [ij, 

Experimentally, QGP search is performed at CERN- 
SPS and RHIC with ultra-relativistic heavy-ion collisions 
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H ES IH El El The experiment of NA-50 collabo- 
ration at CERN-SPS first reported the anomalous J/*5> 
suppression in Pb-Pb collision (158GeV/c per nucleon) 
[l3j . Recently, the RHIC experiments show the various 
signals of QGP, e.g., J/^ suppression 0, enhance- 
ment of stra nge par ticles |16(, jet quenching and highpx 
suppression 0,0] i elliptic flow V2 E3 and so on, in the 
collisions of Au-Au (200GeV/c per nucleon). 

At first, QGP was naively speculated as simple quark- 
gluon gas. Nowadays, there are several indications 
that QGP is not simple perturbative quark-gluon gas in 
quenched lattice QCD. For example, it is pointed out 
that the spatial correlation in tt and a channels remains 
even in QGP phase in lattice QCD E3 ■ The other exam- 
ple is the relation between energy density and pressure, 
which does not satisfy the Stephan-Boltzmann relation 
even above T c (l4|. The calculation of transport coef- 
ficients at finite temperature on quenched lattices also 
shows the strongly correlated gluon plasma [20j. These 
simulations indicate that some of the nonperturbative 
properties may survive in QGP phase. Experimentally, 
in the reports of RHIC UM, GH 03, QGP seems to 
behave as perfect liquid, which strongly interacts, rather 
than dilute gas, from the comparison of the ex peri ments 
with numerical simulations of hydrodynamics 21]. The 
strongly correlated deconfined phase is called strongly 
coupled QGP (sQGP) phase, and is investigated with 
much attention. 

Very recently, some lattice QCD calculations indicate 
an interesting possibility that J/^! and r? c seem to survive 
even above T c [H |H |1 HI . In Ref. jH, the authors 
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calculate correlators of charmonia at finite temperature 
and find the strong spatial correlation between c and c 
even above T c . In Refs. [1313, EH, the authors extract 
spectral functions of charmonia from temporal correla- 
tors at high temperature using the maximally entropy 
method (MEM). Although there are some quantitative 
differences, the peaks corresponding to J/^f and r\ c seem 
to survive even above T c (T c < T < 2T C ) in the cc spec- 
tral function. 

However, all of these calculations may suffer from a 
possible problem that the observed cc state on lattices 
is not a nontrivial charmonium but a trivial cc scatter- 
ing state, because it is difficult to distinguish these two 
states in lattice QCD. One of the reasons of the diffi- 
culty is that an narrow peak does not immediately indi- 
cate a spatially compact (quasi-)bound state. In QGP 
phase, the potential between q and q is considered to be 
the Yukawa potential due to the Debye screening [T3 |. 
Therefore, the binding energy of q and q above T c may 
be small. Then, the bound state of q and q may not be 
spatially compact. In addition, MEM has a relatively 
large error, which sometimes leads to uncertainty for the 
structure of the spectral function. 

In this paper and our previous proceeding |26j . using 
lattice QCD, we aim to clarify whether the cc system 
above T c is a spatially compact (quasi-)bound state or 
a scattering state, which is spatially spread. To distin- 
guish these two states, we investigate spatial boundary- 
condition dependence of the energy of the cc system by 
comparing results on periodic and anti-periodic bound- 
ary conditions. If the cc system is a scattering state, 
there appears an energy difference AE between the two 
boundary conditions as AE ~ 2y/m* + 3tt 2 /L 2 — 2m c 
with the charm quark mass m c on a finite-volume lat- 
tice with L 3 (see Sec.HIjl, If the cc system is a spatially 
compact (quasi-)bound state, the boundary-condition de- 
pendence is expected to be small even in finite volume. 
In Ref. [23, by changing the spatial periodicity of (anti-) 
quarks, the authors actually try to distinguish between a 
scattering state and a spatially compact resonance. 

In this study, we use anisotropic lattice QCD with 
anisotropy a s /a t = 4.0. The reason why we use the 
anisotropic lattice is as follows. At finite temperature 
T, the temporal lattice size is restricted to < t < l/T. 
Then, in calculating the temporal correlator G(t) at high 
temperature, we cannot take the sufficient number of 
points for G(t). By using anisotropic lattice, more data 
points are available for G(t). For the accurate measure- 
ment of G(t), we use such a technical improvement of 
lattice QCD. 

For further technical improvement, we use spatially- 
extended operators with hadron size in the actual lattice 
calculations at high temperature. We are interested in 
the low-lying spectrum in the cc systems at high temper- 
ature in this study, since the ground-state component is 
desired to dominate in the range of < t < l/T. For 
this purpose, we use the spatially-extended operator to 
enhance the ground-state overlap. 



This paper is organized as follows. In Sec. [HJ we 
discuss the method to distinguish a spatially compact 
(quasi-)bound state from a scattering state by chang- 
ing the spatial boundary condition for (anti-)quarks. In 
Sec. IIIII , we briefly explain anisotropic lattice QCD. In 
Sec. II VI we show the method to extract the energy of 
the ground state of cc systems from temporal correlators 
at finite temperature in lattice QCD. Section [V] shows 
lattice QCD results of J/^f and r\ c above T c . Using the 
method discussed in Sec. EI we nn d the survival of J/^ 
and rjc as spatially compact (quasi-)bound states above 
T c (~ 2T C ). In Sec. EU we perform the MEM analysis 
for the cc systems in J '/^S , r\ c and Xci channels above T c 
using the lattice QCD data. Section IVHI is devoted to 
conclusion and outlook. 



II. METHOD TO DISTINGUISH COMPACT 
STATES FROM SCATTERING STATES 

In this section, we explain the method to distinguish 
compact states from scattering states in term of their 
spatial extension. In Sec. Iff Al we discuss and estimate 
the energy shift of a state due to the change of spatial 
boundary condition for the cc systems. In Sec. Ill Bl we 
discuss the correction from a short-range potential. As 
a result, we find that the correction is small compared 
with the energy shift in the cc scattering state above T c . 



A. Boundary condition dependence and energy 
shift for cc scattering states 

For the distinction between compact states and scat- 
tering states, we investigate the cc system on the peri- 
odic boundary condition (PBC) and on the anti-periodic 
boundary condition (APBC), respectively, and examine 
spatial boundary-condition dependence for the cc sys- 
tem. Here, in the PBC and the APBC cases, we im- 
pose periodic and anti-periodic boundary condition for 
(anti-)quarks on a finite-volume lattice, respectively. 

For a compact cc (quasi-)bound state, the wave func- 
tion of the cc system is spatially localized and insensitive 
to spatial boundary conditions in lattice QCD as shown 
in Fig. 1(a), so that the charmonium behaves as a com- 
pact boson and its energy on APBC is almost the same 
as that on PBC [2^. 

For a cc scattering state, the wave function is spatially 
spread and sensitive to spatial boundary conditions as 
shown in Fig. 1(b) and hence there emerges the energy 
difference between PBC and APBC due to the non-zero 
relative momentum of c and c on APBC even in the low- 
est energy state. On PBC, the momentum of a quark or 
an anti-quark is discretized as 

Pk = ^f- (k = 1,2,3, n k eZ) (1) 
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(a) 





Bound state 



Scattering state 



FIG. 1: Schematic figures for boundary-condition dependence 
of a bound state (a) and a scattering state (b). The wave 
function of a bound state is spatially localized, and that of 
scattering state spatially spreads. Therefore, the bound state 
is insensitive to spatial boundary condition, while the scat- 
tering state is sensitive to spatial boundary condition in a 
finite-size box. 



on the finite lattice with the spatial volume L 3 . There- 
fore, on PBC, the minimum momentum is p m i n = 0. On 
APBC, the (anti-)quark momentum is discretized as 



Pk 



(2n k + 1)tt 
L 



(k = 1,2,3, n k G Z). 



(2) 



In this case, the (anti-)quark momentum cannot take zero 
even in the lowest energy state. Then, the minimum 
(anti-)quark momentum is 



\Pn 



V37T 



(3) 



as is depicted in Fig. 2. Thus, there is an energy differ- 
ence between PBC and APBC for the lowest cc scattering 
state. Neglecting the interaction between c and c, the en- 
ergy difference is estimated as 



AE 



2m r 



= 2 v /m2 + 37r 2 / J L 2 -2m c 



(4) 



where m c is charm quark mass. The minimum momen- 
tum of a quark and an anti-quark in the case of scattering 
state is depicted in Fig. 2. 

Table [I] summarizes the mass of the cc compact bound 
state (charmonia) and the energy of the cc scattering 
state both on PBC and APBC. In the case of the cc scat- 
tering state, there emerges the energy difference between 
PBC and APBC. 

Such a method, to distinguish a compact (quasi-) 
bound state from a scattering state by changing bound- 
ary conditions, is actually used in Ref. j2]j. This method 
is essentially based on the fmiteness of lattice volume. 
Note that the finite volume is used in several studies 
for the analyses of a scattering state and/or a compact 
bound state in lattice QCD [2a. I29ll30j| . 






■(tt/L,tt/L,tt/L) 











(-tt/L, -tt/L, -tt/L) 



PBC for cc scattering state 



APBC for cc scattering state 



FIG. 2: Pictorial expression of minimum momentum of the 
(anti-)quark on PBC and APBC in the case of cc scattering 
state in the center of mass flame, c and c are in the finite-size 
with the spatial volume L 3 . On PBC, both particles have zero 
lowest momentum in the lowest state. On the other hand, 
on APBC, quark and anti-quark have an opposite non-zero 
momentum. 



(a) 



S-wave 




FIG. 3: Schematic figure for radial wave functions of cc com- 
pact bound states for the .S'-wave state (a) and the P-wave 
state (b). The S'-wave state can be approximated with a 
Gaussian form e r ' p . In contrast, the wave-function of P- 
wave state should be zero at r = and spatially spreads due 
to the centrifugal potential. 



B. Correction from a short-range potential 

Here, we discuss the possible correction to the energy 
difference of a scattering state between PBC and APBC 
from a short-range potential. In the previous subsection, 
we neglect the interaction between c and c and estimate 
the energy difference between PBC and APBC for the cc 
scattering state. However, in the actual situation above 
T c , a quark and anti-quark interact with each other by 
the Yukawa potential 



V(r) 



-A- 



»(T)r 



(5) 



where m(T) is the temperature-dependent Debye screen- 
ing mass The prefactor A corresponds to the 
Coulomb coefficient at T = and is estimated as A ~ 
0.28 |33|. Therefore, the simple estimation of the en- 
ergy shift AE obtained in Sec. Ill Al may be corrected in 
the presence of the Yukawa potential. In the following, 
considering the Yukawa potential V(r), we estimate the 
energy shift AE of the cc scattering state in the non- 
relativistic quantum mechanics. 
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TABLE I: Summary of the boundary-condition dependence of the compact charmonia and the cc scattering state. The energy 
difference AE is calculated on the spatial lattice size L — 1 .55fm and the charm quark mass m c ~ f .3GeV. 









charmonia 




cc scattering state 


mass 


(energy) on 


PBC 


bound state mass 


M 


~ 2m c 


mass 


(energy) on 


APBC 


bound state mass 


M 




AE e 


= £(APBC) 


- £(PBC) 


~ 




~ Isjml + pi in - 2m c ~ 0.35GeV 



TABLE II: The lowest energy of the cc scattering state above T c on PBC and APBC. Eq c ° and Ej denote the lowest energy 
in the no interaction (free) and the Yukawa-potential cases, respectively (the origin of the energy is shifted by 2m c in this 
non-relativistic estimation). The energy difference from the free case, Eq — Eq cc , is also shown. 





rntrcc 

Eo 




Ea — £jq 


PBC 


OMeV 


-7.8MeV 


-7.8MeV 


APBC 


318.4MeV 


304.7MeV 


-13.7MeV 



Consider two particles, c and c, in the finite box with 
x,y, z £ [0, L], where the boundary condition is periodic 
or anti-periodic, c and c interact each other with the 
Yukawa potential in Eq. JSJ. We estimate the lowest 
energy of the cc scattering state in PBC and APBC cases, 
respectively. 

We use the variational method for the charmonium 
wave function tp(r) (r: the relative coordinate) in this 
estimate. We prepare appropriate basis <fii(r) and expand 
the wave function %j){f) as 



(6) 



where the coefficients {Ci} satisfy J^i \Ci\ 2 = 1, for the 
orthonormal basis 4>i(r). The hamiltonian in the cc sys- 
tem is given by 



H 



W02 02 



d 2 



2/i \dx 2 dy 2 dz 2 



+ V(r), (7) 



where = f v d 3 r <f>\(r)H<f)j(r) and Sij = 

J v d 3 r <fi\{r)(f)j(r). Solving Eq. (fTUfl . the energy of the 
cc scattering state is obtained. 
We take the i-th basis as follows: 



(2 - 5 na ,o)(2 - S ny o)(2 - S nz o) 



L 3 

( 2im x \ ( 2itn v \ ( l-nn 
x cos | — ; — x I cos I — zr^-y I cos I — ; — z 



\ L 



V L 



(11) 



where the normalization factor is set so as to satisfy the 
orthonormal condition, J d 3 r <p\4>j — Sij. In the PBC 
case, we choose the basis as 

ifo :rt= (0,0,0) 

01, h, h-n= (1,0,0), (0,1,0), (0,0,1) 
04, 05, 06 •n = (0,1,1), (1,0,1), (1,1,0) 



where /i — m c /2 is the reduced mass of the cc system. 
The energy of the system is expressed by 



Note that 0i, 02 and 03 degenerate and so on. In the 
APBC case, we choose the basis as 



J v d 3 r ^(f)H^(r) 
J v d 3 r ipi (r)ip(r) 



(8) 



Differentiating Eq. JHJ by Ci and imposing the stationary 
condition for E, i.e., J^P- = 0, we get the equation for E 
and Ci as 



X C i \ I ' f ' r titfHhtf d 3 r 4>\(r)< ! F) 



(9) 



From the condition that {Ci} have non-trivial solutions, 
E is determined as 



det(H ij -ES ij )=0, 



(10) 



& : n 



1 1 1 

2' 2' 2 

3 1 1 

2' 2' 2 
1 3 3 
1 2' 2' 2 



1 3 1 

2' 2' 2 
3 1 3 
2' 2' 2 



1 1 3 

2'2'2 
3 3 1 
2' 2' 2 



In the following, we consider the S'-wave case. Due to the 
spherical symmetry of S'-wave, the wave function -0 is in- 
variant under the replacement of x <-> y, y <-> z, z «-> x. 
Therefore, for instance, the coefficients of 0i,02 and 03 
coincide, and one finds Ci = C2 = C3, C4 = C5 = C*6 
and so on, both in PBC and APBC cases. The conver- 
gence for the number of basis, iVbasis, is checked and the 
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TABLE III: Lattice parameters and related quantities in our anisotropic lattice QCD calculation. 

/3 lattice size a~ L o^ 1 70 "u 3 ut 7f k 
6.10 16* x (14- 26) 2.03GeV 8.12GeV 3.2103 0.8059 0.9901 4.0 0.112 



difference of the energy of lowest state is found to be less 
than lMeV between A^asis = 7 and 13. 

We use A = 0.28 as the Coulomb coefficient in T = 
case. We adopt m(T) = gT as the result of the lowest- 
order perturbative QCD calculation at finite temperature 
[23|. For the QCD coupling constant <?, we use g = 1 as 
a typical value in the infrared region, which corresponds 
to the scale in the typical lattice QCD simulations with 
a — 2Nc ^ q j n ^ n j s ca l cu l a tion, we consider the case 

at T = 300MeV, which is slightly above T c in quenched 
QCD. The mass of charm quark is set by m c = 1.5GeV. 

Table [H] is the summary of the results. i?Q rce denotes 
the energy of the lowest state in free case. Ej denotes 
the lowest-state energy in the Yukawa potential. Note 
that the energy difference \Eq — E f J cc \ is rather small 
as 7.8MeV in PBC case and 13.7MeV in APBC case, 
i.e., about lOMeV. (Even in the extreme case of the 
Coulomb potential, i.e., m(T) = 0, where the potential 
effect is clearly overestimated, the correction is found to 
be 50MeV at most.) 

Thus, the correction from the short-range potential 
about lOMeV is enough small compared with the energy 
shift AE ~ 350MeV between PBC and APBC shown in 
Table [I] Therefore, as far as we consider the cc scattering 
state above T c , short-range interaction between c and c 
can be neglected. 

Here, we note that the parameter set used in this cal- 
culation gives the upper limit of the correction from the 
short-range Yukawa interaction. As the temperature in- 
creases, the Debye screening mass m(T) becomes larger 
and the "Coulomb coefficient" A(T) becomes slightly 
smaller. These tendencies make the interaction smaller 
at high temperature. Thus, using A(T) at T — and 
setting the lowest temperature T ~ T c , which realize 
the QGP phase, we get the upper limit of the correction 
from the short-range interaction. In addition, the use 
of the perturbative estimation m(T) ~ gT also leads an 
"overestimation" for the effect of the Yukawa potential. 
For, the lattice QCD calculation [3lJ shows that the ac- 
tual Debye screening is estimated as m(T) = (2 — 3)gT 
near T c and the qq potential seems to be more screened 
the perturbative result. In fact, the obtained correction 
about lOMeV is reduced in the realistic case. Through 
the above considerations, the correction obtained in the 
calculation can be regarded as the upper limit and the 
actual value of the correction should be smaller than the 
limit about lOMeV. 



III. ANISOTROPIC LATTICE QCD 

In this paper, we adopt anisotropic lattice QCD for 
the study of high-temperature QCD, as was mentioned 
in Sec. Q] We explain anisotropic lattice QCD in the 
following. 

For the gauge field, we adopt the standard plaquette 
action on an anisotropic lattice as [27l l34| 

S ° = 4-^- E ReTr{l-F y ( S )} 

+ 4~ C 1G ? ReTr{l — Pj4(s)}, (12) 

s,z<3 

where denotes the plaquette operator. In the simula- 
tion, we take (3 = 2N c /g 2 = 6.10 and the bare anisotropy 
7c = 3.2103, which lead to the renormalized anisotropy 
as a s /at = 4.0. The scale is set by the Sommer scale as 
Tq 1 = 395MeV. Then, the spatial and the temporal lat- 
tice spacing are evaluated as aj 1 ~ 2.03GeV (i.e., a s ~ 
0.097fm), and a^ 1 ^ 8.12GeV (i.e., a t ~ 0.024fm), re- 
spectively. The adopted lattice sizes are 16 3 x (14 — 26), 
which correspond to the spatial lattice size as L ~ 1.55fm 
and the temperature as (1.11 — 2.07)T C . Here, the critical 
temperature is estimated as T c = (260 — 280)MeV at the 
quenched level [35[- We use 999 gauge configurations, 
which are picked up every 500 sweeps after the thermal- 
ization of 20,000 sweeps. We adopt the jackknife error 
estimate for lattice data. 

For quarks, we use 0(a)-improved Wilson (clover) ac- 
tion on the anisotropic lattice as [27], 03 

S F ee ^2ip(x)K(x,y)ip{y), 
K(x,y) ee 

S x , y K t {{l - li)Ui{x)8 x+l y + (1 +7 4 )02 (a; - 4)^.^} 
-K. J2i(r li)Ui{x)S x+l J + (r + 7i )Uj(x - l)5 x _ ly } 

i 

—k s ce E u iiFii5 Xty — rn s CB E a ijFijdx,y, (13) 
i i<j 

which is anisotropic version of the Fermilab action |36| . 
k s and Kt denote the spatial and the temporal hopping 
parameters, respectively, and r the Wilson parameter, ce 
and Cb are the clover coefficients. The tadpole improve- 
ment is done by the replacement of Ui{x) — > Ui(x)/u s , 
Ui(x) — > U±(x)/ut, where u s and u t are the mean- 
field values of the spatial and the temporal link vari- 
ables, respectively. The parameters k s , Kt,r,CE,CB are 
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FIG. 4: Effective masses of J/* at 1.11T C (a), 1.32T C (b), 1.61T C (c) and 2.07T C (d) in the lattice unit with at = (8.12GeV) _1 . 
The circles and the triangles denote the results on PBC and APBC, respectively. The dashed and dotted lines denote _E(PBC) 
and _E(APBC) obtained from the best-fit analysis, respectively. 



TABLE IV: The energy of the cc system in the J/* channel ( J p = 1") on PBC and APBC at f3 = 6.10 and p = 0.2fm at each 
temperature. The statistical errors are smaller than O.OIGeV. We list also uncorrelated xV-^df and AE = £(APBC)-£(PBC). 



temperature fit ranj 


re S(PBC) [xV^Vdf] 


S(APBC) b?/Nm] 


AE 


I.fLTc 


7-11 


3.05GeV [0.14] 


3.09GeV [0.61] 


0.04GeV 


I.32T C 


8-11 


2.95GeV [0.34] 


2.98GeV [0.33] 


0.03GeV 


1.61T C 


6-9 


2.94GeV [0.10] 


2.98GeV [0.22] 


0.04GeV 


2.07T C 


5-7 


2.91GeV [0.03] 


2.93GeV [0.04] 


0.02GeV 



to be tuned so as to keep the Lorentz symmetry up to 
0(a 2 ). At the tadpole- improved tree-level, this require- 
ment leads to r = a t /a si ce = l/(u s u^), cb = 1/ M s an( i 
the tuned fermionic anisotropy = (utKt) / '(us^s) — 
a s /a t . For the charm quark, we take k = 0.112 with 
1/k = 1/(u s k s ) — 2(~fp + 3r — 4), which corresponds 
to the hopping parameter in the isotropic lattice. The 
bare quark mass mo in spatial lattice unit is expressed 
as too = g(- — 8). We summarize the lattice parameters 
and related quantities in Table 11111 With the present 
lattice QCD setup, the masses of J/'F, r\ c and Xci are 
calculated as toj/^ ~ 3.07GeV, m ric ~ 2.99GeV and 
to Xc1 ~ 3.57GeV at T ~ 0. These values almost the 
same as the experimental ones, mj/^(exp) = 3.10GeV, 
TO r)c (exp) = 2.98GeV and m Xcl (exp) = 3.51GeV. 



IV. TEMPORAL CORRELATORS OF cc 
SYSTEMS AT FINITE TEMPERATURE ON 
ANISOTROPIC LATTICE 

To investigate the low-lying state at high temperature 
from the temporal correlator, it is practically desired to 
use a "good" operator with a large overlap with ground 
state (ground-state overlap) , due to limitation of the tem- 
poral lattice size. To this end, we use a spatially-extended 
operator of the Gaussian type as 

0(t, af) = TV ^ exp c(t, x + y)Tc(t, x) (14) 

with the extension radius p as the hadronic size in the 
Coulomb gauge [13, 13- N is a normalization constant, 
r = 7fe(fc = 1,2,3) and r = 75 correspond to 1~(J/\I / ) 
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TABLE V: The energy of the cc system in the <f] c channel (J p = 0") on PBC and APBC at = 6.10 and p = 0.2fm at each 
temperature. The statistical errors are smaller than O.OlGeV. We list also uncorrelated xV-^df and AE = £(APBC)-£(PBC). 



temperature fit ranj 


;e £(PBC) [xV^Vdf] 


£(APBC) lx z /N DF ] 


AE 


LUTc 


7-11 


3.03GeV [0.04] 


3.02GeV [0.17] 


-O.OlGeV 


1.32T C 


7-11 


2.99GeV [0.78] 


2.98GeV [0.82] 


-O.OlGeV 


i.6ir c 
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FIG. 5: Temperature dependence of the pole mass (en- 
ergy) M J/9 (T) of J/* for (1.11 - 2.07)T C on PBC (circles) 
and APBC (triangles). The circles and the triangles corre- 
spond to _E(PBC) and _E(APBC) of the cc system, respec- 
tively. The energy difference AE = ^(APBC) - £(PBC) ~ 
(0.02 - 0.04)GeV between PBC and APBC are considerably 
smaller than that of cc scattering state, AE ~ 0.35GeV. 



and — (?7 C ) channels, respectively. Note here that this 
form (14) is suitable for S-wave states (see Fig. 3(a)). In 
the actual calculation, the size parameter p is optimally 
chosen in terms of the ground-state overlap. The energy 
of the low-lying state is calculated from the temporal 
correlator, 



G(t) 



(15) 



where the total spatial momentum of the cc system is 
projected to be zero. 

In accordance with the temporal periodicity at finite 
temperature, we define the effective mass m P ff(t) from 
the correlator G(t) by the cosh-type function as [35| 

G(t) cosh[m eS (t)(t~N t /2)} 
G(t + 1) cosh[m ofr (i)(t + 1 -Nt/2)} ' ' 

with the temporal lattice size N t . If the correlator G(t) 
is saturated by the lowest level, it behaves as 



G(t) ~ A cosh 



m (t 



(17) 



where m is the energy of the lowest state. In this case, 
the solution m e g(t) in Eq. (16) coincides with the lowest- 
state energy ttiq in Eq. (17). 



To find the optimal value of the extension radius p, 
we set two criterions. First, an effective mass m c s(t) 
of a temporal correlator G(t) with p should have plateau 
region in terms of t. Second, the ground-state component 
should be large enough in the plateau region of m c g(t). 
For the estimation of the ground state overlap, we define 
9&t(t) as [35| 



g m (t) = A'cosh[m'(t-N t /2)}, 



(18) 



where A' and m' are determined by fitting gat(t) to 
G(t)/G(0) in the plateau region. In general, g&t(0) < 1 
is satisfied, and gat(0) = 1 is achieved if and only if G(t) 
is completely dominated by the ground-state contribu- 
tion in the whole region ofO < t < 1/T. Therefore, 
we use <7fit(0) as an index of the ground-state overlap. 
Namely, if gfit(O) with p\ is larger than that with p2, we 
regard the operator with p\ as a better operator which 
has larger ground-state overlap than that with pi- Un- 
der these criterions, we determine the optimal operator 
in the range of p = (O.lfm) 0.5fm. For p = and O.lfm, 
the effective masses are found to have no plateau region. 
Accordingly, we abandon these extension radii. Next, we 
examine gfit(O) for p = (0.2 — 0.5)fm, and find that the 
optimal size of p is p = (0.2 — 0.3)fm at all tempera- 
tures for cc systems. Actually, the results of p — 0.2fm 
and 0.3fm are almost the same. Hereafter, we show the 
numerical results for p = 0.2fm. 



V. LATTICE QCD RESULTS FOR J/* AND rj c 
CHANNELS ABOVE T c 

We investigate the cc systems above T c both in 
J/^(J = 1~) and r] c (J p = 0~) channels in lattice 
QCD. For each channel, we calculate the temporal corre- 
lator G(t) following Eq. (15), where the total spatial mo- 
mentum is projected to be zero. From G(i), we extract 
the effective mass m e s(t) both on PBC and APBC, and 
examine their spatial boundary-condition dependence. 
From the cosh- type fit for the correlator G(t) in the 
plateau region of effective mass m e g(t), we extract the 
energies, E(PBC) and -E(APBC), of the low-lying cc sys- 
tem on PBC and APBC, respectively. In fact, if the cc 
system is a spatially compact bound state, one can expect 
.E(PBC) ~ £;(APBC), which coincides with the bound 
state mass. On the other hand, if the cc system is a scat- 
tering state, E(PBC) and i?(APBC) are expected to be 
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FIG. 6: Effective masses of rj c at 1.11T C (a), 1.32T e (b), 1.61T C (c) and 2.07T C (d) in the lattice unit with at ~ (8.12GeV) _1 . 
The circles and the triangles denote the results on PBC and APBC, respectively. The dashed and the dotted lines denote 
E'(PBC) and _B(APBC) obtained from the best-fit analysis, respectively. 



the corresponding thresholds of the cc s cattering st ate, 
i.e., £(PBC) ~ 2m c and £(APBC) ~ 2y/m^+p^ in . 

Figure 4 shows effective-mass plots of J/W on PBC 
and APBC for (1.11 - 2.07)T C in lattice QCD. Table llVl 
shows the boundary-condition dependence of the energy 
of the cc system in the J/* channel for (1.11-2.07)T C . ft 
is remarkable that almost no spatial boundary-condition 
dependence is found for the low-lying energy of the cc sys- 
tem, i.e., AE = £;(APBC)-S(PBC) ~ 0. This contrasts 
to the cc scattering case of AE ~ 2^jm? c + 3n 2 /L 2 — 
2m c ~ 0.35GeV for L ~1.55fm and m c ~ 1.3GeV as 
was discussed in Sec. HI This result indicates that J/& 
survives for (1.11 — 2.07)T C as spatially compact (quasi- 
)bound state. Figure 5 shows the temperature depen- 
dence of the pole mass of J/W, Mj/-q,(T). 

Figure 6 shows effective mass plots of r\ c on PBC and 
APBC for (1.11 - 2.07)T C . Table summarizes the cc 
system in the r\ c channel on PBC and APBC at each tem- 
perature. Again, almost no spatial boundary-condition 
dependence is found as AE = £(APBC) - ^(PBC) ~ 0. 
This result indicates that r\ c also survives for (1.11 — 
2.07)T C as spatially compact (quasi-)bound state. Figure 
7 shows the temperature dependence of the pole mass of 
Vc,M Vc (T). 

From Figs. 4 and 6, one finds that m e g(t) on APBC 
is saturated by the low-lying state more rapidly than 
that on PBC, i.e., m^ c (t) > m^ BC (t) for small t 



and m™ c (t) ~ m^ BC (t) for large t. This fact may 
be explained as follows. In the channels of J/vP and r) c> 
there are scattering states in addition to the bound state. 
On APBC, the low- lying scattering states are shifted to 
higher energy region in contrast to the compact bound 
state. Thus, the contributions of scattering states dump 
rapidly on APBC. The (slight) smallness of errorbars on 
APBC compared with those on PBC may be explained 
by the same reason (see Figs. 5 and 7). 

In Fig. 8, we compare M J/9 (T) with M Vc (T) on 
PBC. It is interesting that there occurs the inversion of 
M J/9 (T) and M Vc (T) above 1.3T C . In fact, Mj/^(T) de- 
creases as temperature increases, while M ric (T) is almost 
unchanged (see Fig. 5 and 7). We observe the significant 
reduction of Mj/^(T) of the order of lOOMeV. (We note 
that thermal width broadening also leads to the same 
effect of the pole-mass reduction |23-) In any case, it 
would be interesting to investigate the possible change of 
the J/^f mass above T c . 



VI. MEM ANALYSIS FOR J/*, r] c AND X d 
CHANNELS ABOVE T c 

In this section, we perform the MEM analysis for the 
cc systems in J/'F (J p = 1~), rj c (J p = 0~) and Xd 
(J p = 1 + ) channels above T c using the lattice QCD data. 
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sible P-wave bound state of Xci ■ Actually, we calculate 
the effective mass with Gaussian extended operator with 
p = (0 — 0.5)fm in Xci channel above T c and find that the 
effective mass has no plateau region. This fact indicates 
that we cannot extract the low-lying energy state clearly 
in Xci channel with the operator. 

Instead, we perform the analysis of x r -\ above T c with 
maximally entropy method (MEM) |25| with lo- 

cal interpolating field, Xcl = 0757*0. Using MEM, we 
can extract the spectral function A(u>) from the tempo- 
ral correlator G(t), where the relation between A(w) and 
G(t) at finite temperature T is given by 



FIG. 7: Temperature dependence of the pole mess (energy) 
M Vc (T) of rj c for (1.11 - 2.07)T C on PBC (circles) and APBC 
(triangles). The energy difference AE ~ between PBC and 
APBC are also considerably smaller than that of cc scattering 
state as well as J/'fy case. 



M(T) [GeV] 
3.2 




FIG. 8: Temperature dependence of the pole mass (on PBC) 
of J/* and t) c for (1.11 - 2.07)T C . The squares denote 
Mj/^,(T) and the inverse triangles denote M Vc (T). There 
occurs the level inversion of J/V& and r\ c above 1.3T C . 



Here, the axial vector Xci is a P-wave cc meson and its 
wave function tends to spread due to the centrifugal po- 
tential 1(1 + l)/m c r 2 (r: relative distance between c and 
c, I: orbital angular momentum) compared with S- wave 
states such as J/^S and r\ c (see Fig. 3). According to this 
spread wave function, Xci is expected to be sensitive to 
vanishing of the linear potential and appearance of the 
Debye screening effect above T c . In fact, the dissociation 
temperature of Xd would be lower than that of J /\& and 
j] c . Therefore, we study the P-wave cc meson Xci at finite 
temperature in lattice QCD. 

Here, we note that the difficulty with the extended op- 
erator in axial vector channel. Because Xd is a P-wave 
state, the wave function of Xci is not spherical. More- 
over, its radial wave function should be zero at the origin 
due to the centrifugal potential, as shown in Fig. 3(b). 
Therefore, the spherical and Gaussian type extension of 
the operator, which is suitable for an S*-wave state (see 
Fig. 3(a)), may not be a good choice to investigate a pos- 



G(t) = / duK(t,u)A(u), 
Jo 

K(t,u) = (e-^ + e-( 1 / T -*^)/(l- e -"/ T ). (19) 

We use the Wilson quark action in the calculation of 
lattice QCD data for MEM. Instead of the improvement 
of the action, we use the finer lattice spacing at — a s /4 — 
(20.2GCV)- 1 = 9.75 x 10" 3 fm with = 7.0. The lattice 
size is 20 3 x 46, which corresponds to the lattice volume 
L 3 ~ (0.78fm) 3 and the temperature T = 1.62T C . 

To begin with, we show the MEM results for J/*S> and 
T] c . Figures 9 and 10 show the spectral function of J/* 
and T] c , respectively, with PBC and APBC. Clear low- 
lying peaks corresponding to J/^S and r\ c are observed 
around 3GeV. Here, we note that the appearance of the 
peak structure is highly nontrivial, because the default 
function of MEM is a perturbative one, which does not 
have the peak structure. No difference between PBC 
and APBC is observed for J/^> and r) c , which indicates 
that J/if> and r\ c appear as spatially- localized compact 
bound states. These results of J/\E' and r\ c are consistent 
with those obtained in previous sections. The peaks in 
high-energy region {uj > 5GeV) are considered as lattice 
artifacts |25Ll37l |. 

Now, we investigate the MEM analysis for Xci at fi- 
nite temperature in lattice QCD. Figure 11(a) shows the 
spectral function in Xci(J P = f + ) channel on PBC at 
1.62T C . There is no low-lying peak which corresponds 
to Xci { m xd — 3.5GeV). In fact, in contrast to J/^ 
and T] c , the low-lying structure of the spectral function 
in Xci channel differs from that in J/ty and rj c channels 
at 1.62T C . Therefore, the MEM analysis indicates that 
the dissociation of Xd occurs already at 1.62T C . In the 
high-energy region around 6GeV, we can see a sharp peak 
in the spectral function. Figure 11(b) shows the spectral 
function in Xci channel on APBC at 1.62T C . There is al- 
most no difference between PBC and APBC. We compare 
these results in Fig. 11(c) and confirm that the spectral 
function on PBC almost coincides with that on APBC. 
This BC-independence indicates that the peak around 
6GeV corresponds to a spatially compact (quasi-)bound 
state. The state around 6GeV may be a bound state of 
doubler(s), as was suggested in Refs. [2oll3^ |. 
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FIG. 9: The spectral function A(u>) of the cc state in vector (J p = 1 ) channel on PBC (a) and APBC (b) at 1.62T C extracted 
with MEM from lattice QCD data of the temporal correlator. m(w) denotes the default model of MEM ^ |H El The 
statistical error is denoted by the three solid lines. There is a low-lying peak which corresponds to J/\P (mj/q, — 3.1GeV) even 
above T c . We add (c) as the comparison between PBC and APBC. These results almost coincide. 



VII. CONCLUSION AND OUTLOOK 

We have investigated J/\P and r\ c above T c with 
anisotropic quenched lattice QCD to clarify whether the 
cc systems above T c are compact (quasi-)bound states or 
scattering states. We have adopted the standard Wilson 
gauge action and the 0(a)-improved Wilson quark action 
with renormalized anisotropy a s /a t = 4.0. We have used 
(3 = 6.10 on 16 3 x (14 — 26) lattices, which correspond 
to the spatial lattice volume V = L 3 ~ (1.55fm) 3 and 
T = (1.11 — 2.07)T C . To clarify whether compact char- 
monia survive in the deconfinement phase, we have in- 
vestigated spatial boundary-condition dependence of the 
energy of cc systems above T c . In fact, for low- lying cc 
scattering states, it is expected that there appears a sig- 
nificant energy difference AE = £(APBC) - £(PBC) 
between periodic and anti-periodic boundary conditions 
as AE ~ + 3n 2 /L 2 — 2m c on the finite- volume 

lattice. For enhancement of the ground-state overlap, we 
have used the spatially extended operator with Gaussian 
function. 

As a result, both in J/^ and r\ c channels, we have 
found almost no spatial boundary-condition dependence 
of the energy of the low- lying cc system even on the finite- 
volume lattice for (1.11 — 2.07)T C . These results indicate 
that J 1^1 and rj c survive as spatially compact cc (quasi- 
) bound states for (1.11 — 2.07)T C . Also, the inversion 
of levels of J/^! and rj c above 1.3T C has been seen. In 
fact, we have observed the significant reduction of the 



J/4" pole mass of about lOOMeV, above 1.3T C . Experi- 
mentally, it may be interesting to investigate the possible 
change of the J mass above T c . 

We have also performed the MEM analysis for the 
cc systems in J/\& (J p — 1 — ), r/ c (J p = _ ) and Xci 
(J p = 1 + ) channels above T c using the lattice QCD 
data. For this analysis, we have adopted lattice QCD 
at 1.62T C with the Wilson quark action with j3 = 7.0 
and a s /at = 4.0. For the 5-wave channel, we have ob- 
tained the same results on the survival of J and r\ c as 
the spatially- localized compact bound state above T c . In 
contrast to the J and r\ c cases, the spectral function 
in the P-wave channel has no low-lying peak structure 
corresponding to Xci around 3.5GeV. This fact indicates 
that Xci already dissociates at 1.62T C . The further anal- 
ysis of charmonia at high temperature with MEM will be 
reported in Ref. [38l| . 

Through the MEM analysis, we have observed the spa- 
tially localized bound state in high energy region, which 
is not affected by the spatial boundary condition at all. 
The bound state appearing in high-energy region on the 
lattice would be the bound state of doubler(s), which is 
unphysical. 

Our study of charmonia indicates that the S-wave 
mesons J/ty and r\ c survive even above T c (~ 2T C ) as com- 
pact bound states. 

As a successive work, we are performing the same anal- 
yses for other charmed mesons, e.g., D mesons. The 
narrowness of the decay width of charmonia strongly de- 
pends on whether the decay channel of DD opens or not. 
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FIG. 10: The spectral function A(u>) of the cc state in pseudo scalar (J p = _ ) channel. The notations are the same as those 
in the J/$ case (Fig. 9). There is a low-lying peak which corresponds to rj c (m Vc ~ 3.0GeV). Because the results in PBC and 
APBC cases almost completely coincide, the dotted line is blinded by the solid line in Fig. 10(c). 



If the mass of the D meson is shifted at high temperature, 
a drastic change of decay width of some charmonia pos- 
sibly occurs, which may be interesting both theoretically 
and experimentally 39]. 

The subject of hadrons in high density system is fas- 
cinating. For example, CERES collaboration presented 
interesting events on the lepton pair production in low- 
mass region from p-Be and p- Au ( 450GeV/c) and from 
S-Au (200GeV/c per nucleon) |4(J. It is reported that 
the lepton pairs from S-Au in the invariant mass range 
0.2GeV/c 2 < m < 1.5GeV/c 2 are largely enhanced com- 
pared with the hadronic contributions, while those from 
p-Be and p-Au are not enhanced. The enhancement is 
considered to be the signal of the change of hadronic na- 
ture of vector mesons in finite density system |40L l4ll |42| . 
Namely, the mass reduction and/or the width broaden- 
ings of vector mesons at finite density can lead to the 
enhancement of low-energy lepton pairs. At the exper- 
imental facility J-PARC in Japan, which is planned to 
run in 2008, the experiments of high density system are 
expected to take place. For these experiments, the the- 



oretical study of finite density system is required. Then 
the study at finite density in lattice QCD is challenging 
and worth trying. 

The survival of J/\I> above T c may change the scenario 
of J suppression. These analyses give us the further 
knowledge of QCD at finite temperature. 
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